Collective luminescence and phonon-induced processes in double quantum dots 
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We study the evolution of a quantum state of a double quantum dot system interacting with 
the electromagnetic environment and with the lattice modes, in the presence of a coupling between 
the two dots. We propose a unified approach to the simulation of the system evolution under joint 
impact of the two reservoirs. We discuss the sub- and superradiant radiative decay of the system, 
the phonon-induced decay of entanglement between the dots, and the transfer of excitation between 
them. 
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I. INTRODUCTION 

Semiconductor structures composed of two closely 
spaced quantum dots (QDs) have attracted enormous at- 
tention in recent years. Pairs of vertically stacked QDs 
with spatial separation down to single nanometers can 
be obtained in a two-layer self-assembled growth pro- 
cess, where the strain distribution favors QD nucleation 
in the second layer on top of the QDs formed in the first 
layer [l], . The development of the manufacturing tech- 
nologies [3, H| has made it possible to achieve structures 
built of nearly identical dots, with the splitting of ground 
state transition energies down to single meV (very likely 
interaction- limited) [f| @] . The state space of such a dou- 
ble quantum dot (DQD) is obviously richer than that of 
a sing le QD @,| and allows, e.g., for entanglement be- 
tween the dots. Also the recombination and relaxation 
processes in DQDs show many features which cannot ap- 
pear in individual QDs. The quantum coherence of car- 
rier states in DQDs are affected by the interference and 
collective effects that appear in the interaction of such 
systems with their radiative environment (electromag- 
netic vacuum) and with the surrounding crystal lattice 
(phonons). 

Optical properties of DQDs may be strongly modified 
due to collective interaction of sufficiently closely spaced 
QDs with the electromagnetic (EM) field. These col- 
lective effects have been extensively studied for atomic 
systems d where they manifest themselves by superra- 
diant emission [ic| . A signature of superradiant behavior 
was also observed in ensembles of QDs [ll| . On the other 
hand, the collective interaction leads to the appearance of 
subradiant states which are decoupled from the environ- 
ment and, therefore, do not undergo decoherence. It has 
been proposed to use these states for noiseless encoding 
of quantum information [l2l [l3j . 

Another area where new features appear for DQD sys- 
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terns is the decoherence due to carrier-phonon coupling. 
Dephasing of carrier states in QDs due to carrier-phonon 
dynamics has been observed experimentally as a decay 
of nonlinear optical response in a four-wave mixing ex- 
periment with ultrashort pulses [3, UH . A characteristic 
feature of the phonon-induced dephasing in QDs is that 
it is only partial, i.e., after a few picoseconds of carrier- 
phonon dynamics, the degree of coherence (i.e., the val- 
ues of the off-diagonal elements of the density matrix) 
reaches a certain finite level, depending on the system 



geometry and temperature [lei, ]. In DQDs and regu- 
lar QD arrays, the degree of dephasing may be reduced 
by encoding the logical qubit values into many-exciton 
states over a QD array [l8|. Phonon-induced dephasing 
is also detrimental to entanglement in DQDs and larger 
QD arrays. The impact of the partial dephasing on the 
entanglement is very strong since the latter is more prone 
to dephasing than local coherence and may be completely 
destroyed even though the decoherence is only partial 
[l9l I20I [2lj . Moreover, because of the delicate aspect of 
inter-subsystem coherences involved, the decay of entan- 
glement due to dephasing strongly depends on the na- 
ture of environmental interaction (the same vs. different 
reservoirs for the two subsystems) (22j . 

Various types of coupling between the dots change the 
properties of the system even further. Theoretical calcu- 
lations show that for closely spaced dots, tunnel coupling 
(wave function overlap) between the dots should strongly 
affect their electronic structure [23|, [24], El . Optical spec- 
tra of such structures indeed show clear manifestations 
of electronic coupling @, El, E3, El El • 

On the other hand, wave function overlap is not the 
only mechanism of interaction between the QDs. In fact, 
for QD separations of about 10 nm the energetically low- 
est states (in the absence of external fields) correspond 
to spatially direct excitons localized in individual QDs 
[25| . Such states are still bound by Coulomb interaction. 
While the static ( "direct" ) dipole coupling preserves the 
occupations of the individual QDs, the Forster interac- 
tion via interband dipole moments [3(J l3lll (first intro- 
duced in the context of molecular systems [32|, [33|) makes 
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it possible to transfer the exciton occupation between the 
dots. In contrast to the tunnel coupling, which is analo- 
gous to chemical bonding between atoms and, therefore, 
turns the two dots into one quantum system (a quan- 
tum dot molecule), the dipolc couplings are rather like 
van der Waals forces between separate entities. There- 
fore, QDs coupled by this kind of interaction are much 
closer to the general paradigm of well-defined, separate 
qubits [H| on which, however, collective quantum opera- 
tions (multi-qubit gates) can be performed, based on the 
interactions between the systems. 

In a closed system the (usually very weak) Forster in- 
teraction has considerable effects only very close to res- 
onance [30L I3II]. However, carrier-phonon coupling pro- 
vides the necessary dissipation channel which opens a 
possibility of excitation transfer driven by the Forster in- 
teraction even if the energy mismatch between the dots 
is much larger than the interaction energy. Phonon- 
assisted excitation transfer between the quantum states 
of a molecule was in fact observed in many experiments 
BUI, S S3, Hill- In closely stacked dots this excita- 
tion transfer process is mostly due to phonon-assisted 
tunneling of carriers. However, for larger separations 
tunneling is exponentially suppressed. Indeed, phonon- 
assisted transitions involving tunneling are very ineffi- 
cient for a 10 nm separation even though a small energy 
splitting matches the acoustic phonon energies [39| . In 
such cases, the transfer is most likely due to the Forster 
coupling. 

This paper presents a theoretical review of a few phe- 
nomena that appear in double-dot systems. The existing 
results, which we originally obtained using various theo- 
retical approaches, are presented within a unified formal- 
ism. Moreover, some examples of the interplay between 
the effects due to the phonon and photon reservoirs are 
discussed. In Sec. |H] we define the model and our ap- 
proach to the simulation of its evolution. In Sec. IIIII 
we present a theory of excitation transfer between cou- 
pled quantum dots due to interband dipole coupling and 
carrier-phonon interaction [4(j. Next, in Sec. II VI we sum- 
marize our studies on collective effects in spontaneous 
emission [4l|. We show that superradiant-like behavior 
may appear in the optical response of sufficiently strongly 
coupled pairs of quantum dots even if the two dots are 
not identical. In addition, we extend the existing result 
by discussing the effect of phonon-induced transitions on 
the collective optical properties. In Sec.[Vl we discuss the 
decay of entanglement due to phonon-induced dephasing. 
Again, we extend the existing study of phonon-related 
effects on picosecond time scales [42| by discussing the 
long-time decay due to spontaneous emission and the 
mutual impact of phonon-related and radiation-related 
effects. The final Sec. I VII contains concluding remarks 
and some outlook for a possible further development of 
the theory. 



II. THE SYSTEM 

The system under study is composed of two coupled 
QDs with transition energies t\ and £2, interacting with 
their phonon and photon (radiative) environments. We 
restrict the discussion to the ground states of excitons in 
each dot and assume that the spin polarizations of the 
excitons are fixed. As the exciton dissociation energy in 
absence of external electric fields is rather large (several 
to a few tens of meV), we consider only spatially direct 
exciton states, i.e., such that the electron-hole pairs re- 
side in one and the same dot. 

Thus, the model includes four basis states |mn), 
to, n = 0, 1, where to and n denote the number of excitons 
in the first and second dot, respectively. We will use ei- 
ther the explicit tensor product notation or a contracted 
one both for the states of the two dots (\mn) = |to)<8) |n)) 
and for the operators (\mn)(m'n'\ = |to)(to'| <g) \n)(n'\). 

We will describe the evolution in a "rotating basis" 
defined by the unitary transformation 

jj _ e iBt(|01>(01|+|10>(10|+2|H){U|)/R 

where E = (ei + e-z)/2. The Hamiltonian is then 

H = Huqd + H p h + i/ ra d + -ffc-ph + #c-rad- 

The first term describes exciton states in the DQD 
structure. 

ff DQ D = A(|1)(1|(8I-I(8|1)(1|) 

+F(|l)<0|®|0)<l|+H.c), (1) 

where A = (ei — e 2 )/2, 1 is the unit operator, and V is the 
coupling between the dots, which can be assumed real. 
Electron and hole wave functions will be modelled by 
identical anisotropic Gaussians with identical extensions 
I in the xy plane and l z along z for both particles, 

/ 1 x 2 + y 2 lz 2 \ 
^(r)-exp -_). 

The coupling V may originate either from the Coulomb 
(Forster) interaction or from tunnel coupling. In the for- 
mer case it is related to the distance between the dots 
and to band structure parameters [13, El[ , 

e 2 \a\ 2 

Here e is the electron charge, eo and e r are the vacuum 
permittivity and the dielectric constant of the semicon- 
ductor, respectively, and 

hP cv h 
a = w — . 

TOo-E-g ^/2m e E g 

where P cv is the interband marix clement of the momen- 
tum operator, too and m e are the free and effective elec- 
tron masses, and E g is the band gap. The function f(x) 
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accounts for the correction to the point dipole approxi- 
mation due to the finite size of the dots (in fact, compa- 
rable to D) 0, S3 and is given by 



u(t) - x 2 t 2 
u 5 / 2 {t) 



dt{\ - t 2 ) 



■ exp 



2u(t) 



where u(t) — 1 — t 2 + (l z /l) 2 t 2 . This correction reduces 
the coupling for D < I and removes the 1 /D 3 singularity 
at D — ► 0, while it does not affect the coupling for D 3> I 
since f(x) — > 1 as x — > oo. We will neglect the possible 
biexciton shift (coupling between static, intraband dipole 
moments). 

The phonon modes are described by the free phonon 
Hamiltonian 



H, 



ph 



where bk,b\, are bosonic operators of the phonon modes 
and LUk are the corresponding frequencies. Interaction of 
carriers confined in the DQD with phonons is modelled 
by the Hamiltonian 

ff c _ ph = (|l)<l|®I)£/«(4 + 6_ fe ) (3) 
k 

+(I®|lXl|)£/f (b{ + b- k ), 

k 

(12) 

where ' are system-reservoir coupling constants. For 
Gaussian wave functions, the coupling constants for the 
deformation potential coupling between confined charges 

(1 2) 

and longitudinal phonon modes have the form = 
fke ±lkzD ^ 2 , where D is the distance between the subsys- 
tems and 



fk = (o"c — fh)l 



2gvc\ 



exp 



Here v is the normalization volume, k±/ z are momentum 
components in the xy plane and along the z axis, er e /b 
are deformation potential constants for electrons/holes, 
c\ is the speed of longitudinal sound, and g is the crystal 
density. 

The third component in our modeling is the radiative 
reservoir (modes of the electromagnetic field), described 
by the Hamiltonian 



H, 



rad 



Hwkcl x c k,\, 



where Ck,\, c k A are photon creation and annihilation op- 
erators and Wk are the corresponding frequencies (A de- 
notes polarizations). The QDs are separated by a dis- 
tance much smaller than the relevant photon wavelength 
A = 2irhc/E, where E = (ej — €2) [2, so that the spa- 
tial dependence of the EM field may be neglected (the 
Dicke limit). The Hamiltonian describing the interaction 



of carriers with the EM modes in the dipole and rotating 
wave approximations is 



(4) 



with 



and 



= |0)(l|(g)I + I(g)|0)(l| 



g k \ = id ■ a\(k) 



2eoe r v ' 



where d is the interband dipole moment (d = ea) and 
e\(k) is the unit polarization vector of the photon mode 
with the wave vector k and polarization A. For wide- 
gap semiconductors with E ~ 1 eV, zero-temperature 
approximation may be used for the radiation reservoir at 
any reasonable temperature. 

In certain limiting cases, analytical formulas for the 
evolution of the DQD system may be found. For uncou- 
pled dots (V = 0) interacting only with lattice modes 
(phonons), an exact solution is available If only 

the radiative decay is included, a solution in the Markov 
limit can be obtained [4JJ . Here, we will use a description 
which allows one to deal with the simultaneous action of 
both these environments. We describe the evolution of 
the reduced density matrix of the DQD system in the in- 
teraction picture with respect to -/?dqd by the equation 

P = A-ad[p] + £ P h[p\- 

Here the first term describes the effect of the radiative 
decoherence in the Markovian limit in terms of the Lind- 
blad dissipator 



£ ra d [p] = r 



rad 



1 



£_(t)p£+(t) --{£+(£)£_(*), p}_ 



where £_(i) = E^(t) = e !H DQDt/s s _ e -fflDQD*/s and 

2 



Trad — 



E 3 \d\ W e r 



3ne c 3 h 4 



is the spontaneous decay rate for a single dot. The second 
term accounts for the interaction with the non-Markovian 
phonon reservoir. We use the time-convolutionless equa- 
tion 

£ P h[p] = - dr Tr ph [i?c-ph(*), [-ffc-ph(r), p(t) eg p ph ]) , 
Jo 

(5) 

where H c ^ ph {t) = e %HD » Bt ' H H c -p} l e~' tHDqiDi/h , Pph is the 
phonon density matrix at the thermal equilibrium, and 
Tr p h denotes partial trace with respect to phonon degrees 
of freedom. 

The above equation of motion for the reduced density 
matrix strictly reproduces the results in the limiting cases 
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FIG. 1: The occupation of the higher energy QD as a function 
of time at T — 4 K: (a) D — 8 nm and A as shown; (b) 
A = 3 meV and D as shown. 



mentioned above. Moreover, for the case of a DQD cou- 
pled to phonons with non- vanishing coupling V, it yields 
results reasonably close to those obtained by a much more 
complex correlation expansion technique [40(. 

In numerical simulations, we take the parameters 
corresponding to a self-assembled InAs/GaAs system: 
(j c — Oh = 9 eV, p = 5350 kg/m 3 , c\ = 5150 m/s, 
the wave function parameters I = 4.5 nm, l z = 1 nm, 
and the radiative recombination time (for a single dot) 
l/r rad = 400 ps. 



III. PHONON-ASSISTED EXCITATION 
TRANSFER 

In this section we discuss the evolution of a DQD sys- 
tem coupled only to its phonon reservoir. The system is 
then described by the Hamiltonian H = i?DQD + #ph + 
H c - P h, [Eqs. |T]) and ©]. We will assume that there is 
one exciton in the system, initially localized in one of the 
dots. We will see that the interplay of the coupling and 
phonon-assisted dissipation leads to irreversible excita- 
tion transfer between the dots, which is due to a weak 
Coulomb (Forster) coupling between them. 

In general, if the initial state corresponds to the ex- 
citon located in one of the dots (state 1 01) or 1 10} ) the 
evolution is a combination of conservative (unitary) oscil- 
lations due to the perturbation induced by the coupling 
V and a dissipative, irreversible transition towards the 
lower eigenstate of Hdqd [S EH induced by the cou- 
pling to the phonon continuum. However, in the case 
of a weak coupling the oscillations are very small and 
the occupation of the higher-energy dot follows a nearly 
exponential decay, as shown in Fig. [TJ In these compu- 
tations, the coupling energy V is calculated according to 
Eq. An interesting feature visible in Fig. [1] is that 
the rate of the excitation transfer is non-monotonic both 
in the energy mismatch and in the separation between 
the dots. 

The nearly exponential decay curve suggests that the 
process can be described in the Markovian approxima- 
tion. This is possible, since the typical time scale of 
the decay process is long compared to the transition fre- 
quency V/h between the two single-exciton states. We 



FIG. 2: The rate of phonon-assisted excitation transfer as a 
function of the energy mismatch for a few values of the QD 
separation D (a) and as a function of D for a few values of A 
(b) at T = 4 K. 

have shown [40| that the formal long-time limit of Eq. (J5]) 
in the rotating wave approximation leads to optical Bloch 
equations describing the system dynamics in the inter- 
action picture and in the rotating frame related to the 
eigenstates of -£/dqd- 

In general, upon transforming back to the original ba- 
sis 1 01 ) , 1 10) and to the Schrodinger picture one obtains 
a complex evolution. However, in the case of |V| <C |A|, 
which is of particular practical importance, the eigen- 
states of i?DQD are very close to 1 10} and |01). More- 
over, in this limit the energy difference between the eigen- 
states, Ml = 2\ / V 2 + A 2 is nearly equal to 2A. In this 
limit, one obtains an exponential excitation transfer with 
the rate 

r = 47r^ [#(fi) + R(-Q)] , 

with the spectral density 

R(w) = | |n B (w) + sin 2 ^ \f k \6(\u\ - w k ). 

k 

Thus, the Markovian equations are particularly useful in 
the limit of weak coupling, where the Markovian dephas- 
ing rate may be identified with the rate of irreversible 
excitation transfer between the dots. 

The rate for the phonon-assisted process is governed, 
on the one hand, by the amplitude of the Forster coupling 
which decreases roughly as 1/D 3 . On the other hand, it 
is strongly influenced by the structure of R(£l) . The lat- 
ter oscillates as a function of both and D [4fj| due to the 
interplay between the wavelength of the emitted phonon 
and the QD separation in the molecule (phonons are pref- 
erentially emitted along the strongest confinement limit, 
i.e., along the DQD axis). In particular, R(fl) has a pro- 
nounced minimum whenever Q is a multiple of 2-kHc\/ D 
which explains the oscillating dependence on both A and 
D (see Fig. [2). 

The transfer rate for A = 4 meV, D = 4.5 nm is 
r = 2.3 ns -1 , which is about twice lower than the value 
of r = 5.25 ns -1 deduced from fitting to the photon 
correlation data 0. It should be noted, however, that 
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FIG. 3: The exciton occupation for sub- (a) and superradiant 
(b) states for A = 1 meV. The inset in (b) shows the values 
of the occupation decay rates. Line definitions in (a) refer to 
both plots. 



our modelling is based on certain choices of parameters 
that cannot be uniquely determined. First of all, the 
magnitude of the Forster coupling for a given inter-dot 
distance can only be roughly estimated, since the value of 
the parameter a in an inhomogeneous, strained structure 
is not exactly known. Moreover, the carrier-phonon cou- 
pling constants in the relevant energy (or wave vector) 
range are strongly geometry dependent. For instance, by 
changing the localization widths to l c — lh — 4.0 nm and 
l z = 0.8 nm one gets a considerably increased value of 
r = 3.65 ns" 1 . 



IV. COLLECTIVE RADIATIVE PROPERTIES 
OF DOUBLE QUANTUM DOTS 

In this section, we study the effect of the collective 
coupling between the two interacting QDs and their elec- 
tromagnetic environment on the radiative relaxation of 
carriers. We begin with a model of a DQD interacting 
only with the EM field. Then, we discuss the effect of 
the additional phonon-induced relaxation. 

In the case of purely electromagnetic environment, the 
Hamiltonian of the system is H = ifDQD+#rad+-#c-rad- 
This Hamiltonian conserves the number of excitations 
(excitons plus photons). In this case, the solution ob- 
tained using the equation ^ is equivalent to that derived 
using the Wigner-Weisskopf approach [il[ . 

In Fig. [3] we show the evolution of the average number 
of excitons in the DQD for two initial states 

|tf ± (0)> = (|01>±|10>)/V2, 

for a DQD with a realistic value of the energy mismatch 
A = 1 meV. For V <C A, both states show simple expo- 
nential decay with the rate r ra d- This is understandable, 
since two different dots emit radiation into different fre- 
quency sectors of the reservoir and no collective effects 
should be expected. The situation changes in the oppo- 
site limit, V ^> A. Now, one of the states becomes stable 
(subradiant), while the other (superradiant) state decays 
exponentially with a twice larger rate. It should be noted 
that the energy mismatch is much larger than the emis- 
sion line width, so that the subradiance and superradi- 
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FIG. 4: The exciton occupation for sub- and superradiant 
states in the presence of phonon-induced relaxation for A = 
1 meV, d = 8 nm, T = 4 K. (a) Subradiant state for V 

1 meV; (b) superradiant state, V = 1 meV. Red solid lines 
show the evolution without phonons and the green dashed 
lines with phonons. 



ance effect is due to the coupling between the dots. In the 
intermediate range of parameters, the decay is not expo- 
nential. It can be shown [4ll | that for the superradiant 
initial state, the number of excitons evolves as 



n(t) = \c(t)f 



where 



and 



sin ip 



1 

71 



+ cos 2 ((^ + V4)e 2RcA+t , 

1/2 



VA 2 + V< 



A± = 



Trad 



± V-a 2 + W + r rad /2) 5 



The values of the two exponents X± are shown in the 
inset of Fig. [3(b). 

The subradiance and superradiance effects discussed 
above appear in spite of the large energy mismatch be- 
cause the single-exciton eigenstates of the DQD Hamilto- 
nian are superpositions of the basis states 1 01} and 1 10} 
which are partly of sub- or superradiant character (for 
the —1 and +1 relative phase, respectively). It is there- 
fore clear that the stability of these states is essential for 
the effect. 

When carrier-phonon interaction is included, these 
states are no longer stable due to excitation transfer pro- 
cesses discussed in Sec. IIII1 Now, the exciton popula- 
tion decay will depend on the interplay of radiative and 
phonon-related effects. This is illustrated in Fig. [4] where 
we plot the exciton number as a function of time for the 
same two initial states as previously, but now we include 
the coupling to phonons. Contrary to the purely optical 
case, the kinetics now depends on the sign of the cou- 
pling. For V < 0, the lowest single-exciton state has a 
superradiant character. Due to phonon-induced relax- 
ation, the system undergoes transition to this state and 
the subradiance effect is partly destroyed [see Fig. UK a)]. 
An opposite situation takes place for V > [Fig. ffljb)]. 
Now, the phonon-induced relaxation consists in a tran- 
sition to the subradiant state. As a result, the radiative 
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recombination is slowed down by the additional phonon- 
related decoherence. 

Let us note that the superradiant state is particularly 
relevant for optical experiments since such a bright com- 
bination of single exciton states is excited by ultrafast 
optical pulses from the ground state. Thus, the curves 
plotted in Figs.^b) andrjjb) directly correspond to the 
decay of population after an optical excitation. 



ENTANGLEMENT DECAY 




FIG. 5: Evolution of entanglement for uncoupled dots (V = 0) 
at various temperatures. 



In this section we discuss the evolution of entanglement 
between the states of two QDs under dephasing caused 
by carrier-phonon and carrier-photon interactions. In 
our earlier work [42| we studied two dots in the absence 
of excitation transfer coupling, undergoing only phonon- 
induced pure dephasing. Since that problem was repre- 
sented by an independent boson model the evolution of 
the open system could be found in a closed analytical 
form. We showed that phonon-induced pure dephasing, 
in spite of its only partial character, can lead to com- 
plete decay of entanglement after a finite time. As we 
pointed out, this happens for certain initial maximally 
entangled states in which all the coherences are present, 
that is, all four basis states are involved in the superpo- 
sition. We studied also the dependence of the entangle- 
ment decay on the distance between the dots. We showed 
that the degree of dephasing-induced disentanglement in- 
creases with growing distance between the dots and that 
non-zero distance is a necessary condition for complete 
disentanglement. As expected, the degree of dephasing 
and, therefore, disentanglement increases as the temper- 
ature grows. The decay of entanglement due to coupling 
with the electromagnetic field was also studied [20l |. 

Here, we extend our earlier model [HJ by including 
the radiative decay of the entangled excitons (coupling to 
the photon reservoir) and excitation transfer interaction 
(V 7^ 0) between the dots. Since the extended model be- 
comes more complicated and presents many new features, 
we focus on these new elements and restrict the discus- 
sion to one initial state and fixed distance D — 8 nm 
between the dots. We assume also that the biexciton 
shift is absent. The dots are assumed different, with the 
energy mismatch of A = 1 mcV. In the presence of spon- 
taneous emission and inter-dot coupling, the evolution 
equation cannot be solved analytically and we use the 
numerical approach described in Sec. HU We study the 
evolution of a maximally entangled "singlet" state 



m = o)) = 



|01> - |10> 
V2 ' 



As a measure of entanglement, we use the Wootters con- 
currence [H,[I3- 

The evolution of the entanglement between the dots 
is shown in Fig. [5] At short (picosecond) time scales, 
phonon-induced dephasing leads to a drop of concur- 
rence, ending with a temperature-dependent plateau 




FIG. 6: Evolution of entanglement for coupled dots (V = 
±1 meV, as shown) compared to the uncoupled case (V = 0) 
at T = 40 K. 



level. In the absence of radiative recombination, the en- 
tanglement would remain constant after this initial de- 
phasing stage. However, in the presence of carrier-photon 
coupling the exciton lifetime becomes finite, which leads 
to an exponential population decay and, in consequence, 
to decay of entanglement on the time scales ~ l/r rac j. 

The situation becomes much more complicated (and 
more interesting) if the dots are coupled by a transfer- 
type interaction (that is, V ^ 0). The evolution in this 
case is shown in Fig. A few effects can be seen. The 
most striking feature are the oscillations of concurrence. 
Since the coupling is comparable to the energy mismatch, 
the system performs rotations in the single-exciton sub- 
space, coming close to the separable states 1 01) or 1 10) 
every half-period. These oscillations are damped on a 
time scale of tens of picoseconds, as the excitation is dis- 
sipatively transferred to the lower-energy eigenstate of 
Hbqd by a process discussed above in Sec. IIIII Depend- 
ing on the sign of the interaction, this eigenstate (which 
can still be entangled) can have either subradiant or su- 
perradiant character (for V > and V < 0, respectively). 
This is visible as the difference in the entanglement decay 
rates between the two cases for long times. 



VI. CONCLUSIONS AND OUTLOOK 

In this paper, we have reviewed and extended our re- 
cent results on the interaction between systems composed 
of two quantum dots and their environment. We have 
shown that collective interaction between the dots and 
the surrounding radiation reservoir leads to sub- and su- 
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perradiance effects. In the presence of a coupling be- 
tween the dots, these collective effects become more sta- 
ble against the differences between the two transition 
energies. We have also seen that carrier-phonon cou- 
pling leads to new dephasing effects in such double-dot 
systems. In particular, phonon-induced decoherence de- 
stroys entanglement between the two dots. This destruc- 
tion of entanglement is much stronger than local dephas- 
ing. We have discussed the joint effect of phonon-induced 
dephasing and inter-dot coupling. In this case, the most 
important feature of the system dynamics is the irre- 
versible excitation transfer between the dots. 

A new area of interesting and important phenomena 
emerges if both the spontaneous emission and phonon- 
induced dephasing are simultaneously taken into account. 
Most of this field remains to be exploited. As an example 
of physical effects that may appear due to joint impact of 
the two reservoirs we have discussed the phonon-induced 
modification of collective emission due to transitions be- 
tween subradiant and superradiant states and the inter- 
play of the two dephasing channels in the decay of en- 



tanglement between excitons confined in the two dots. 

Obviously, the phenomena discussed here are only a 
fraction of all the physical effects that may be present 
in a system of two non-identical quantum dots interact- 
ing with two different reservoirs (photons and phonons) 
in the presence of coupling between the dots. In partic- 
ular, the evolution of double quantum dot systems un- 
dergoing simultaneous dephasing via both radiative and 
phonon-related channel has been studied only in a very 
limited scope. Investigation of these problems will cer- 
tainly bring new knowledge not only on the properties of 
quantum coherence in this specific semiconductor system, 
but also on general properties of open quantum systems. 
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